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Preface 


(  This  study  was  undertaken  for  the  Control  System 
Development  Branch,  Air  Force  Flight  Dynamics  Laboratory 
(Flight  Control  Division).  The  first  objective  of  the  effort 
was  to  ascertain  whether  the  existing  linear  mathematical 
model  of  the  AC130  was  accurate.  I  thought  it  best  to  begin 
with  the  basic  equations  of  motion  and  develop  a  parallel 
model,  the  characteristics  of  which  could  be  compared  to 
those  of  the  existing  model. 

Using  the  model  developed  in  this  study,  the  second 
purpose  was  to  compile  date  on  selected  stability  and  control 
parameters  for  the  AC 130 A  aircraft  to  predict  general  trends 
and  values  for  these  parameters  throughout  an  extensive 
flight  envelope. 

I  wish  to  express  my  gratitude  to  the  following 
personnel:  to  Lt.  Col.  James  Thompson,  for  sponsoring  this 
study;  to  Capt.  Kenneth  Bassett,  for  supporting  the  project 
and  guiding  my  efforts  to  produce  useful  results;  to  Capt. 

Len  Kruczynski,  for  supplying  the  basic  computer  routines 
used  in  the  computations;  and  to  my  wife,  for  her  unending 
patience  with  the  author. 


Robert  G.  Lorenz 
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Notation 


Symbol 

• 

•  Description 

Units 

A 

Altitude 

ft 

a 

Speed  of  sound 

ft/sec 

A.C. 

Aerodynamic  center,  %  MAC 

AR 

Aspect  ratio 

b 

Span 

ft 

C.G. 

Center  of  gravity,  %  MAC 

«# 

Section  .chord 

ft 

c 

j  • 

Mean  aerodynamic  chord  (MAC) 

ft 

* 

9 

D 

Drag 

lbf 

<1 

Diameter  of  propeller 

ft 

F 

Flap  position  in  % 

8 

Acceleration  due  to  gravity 

ft/sec2 

Scaler  components  of  angular 

2 

a  y  s 

momentum  of  spinning  rotors 

slug  ft  /sec 

I 

Total  moment  of  inertia  of 
propeller 

slug  ft2 

x« 

Moment  of  inertia 

slug  ft2 

iy 

Moment  of  inertia 

slug  ft2 

*. 

Moment  of  inertia 

slug  ft2 

Ixz 

Product  of  inertia 

slug  ft2 

*xy 

Product  of  inertia 

slug  ft2 

Xyz 

Product  of  inertia 

slug  ft2 

i 

Angle  of  incidence 

deg 

K 

Induced  drag  coefficient,  1/weAR 

L 

Lift 

lbf 

•  • 
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Roiling  moment  about  mass  center 


lbf  ft 


Longitudinal  distance  between  center 
of  gravity  and  of  aerodynamic  center 
hbrizontal  tail  (assumed  to  be  at  .25c) 

If  Longitudinal  distance  between  center 

*  of  gravity  and  of  aerodynamic  center 

vertical  tail  (assumed  to  be  at  .25c) 

M  Mach  number 

MAC  Mean  aerodynamic  chord,  c 

(ft  Pitching  moment  about  mass  center 

m  Mass 

(t  Yawing  moment  about  mass  center 

P,Q,R  Scaler  components  of  angular  velocity 

p,q,r  Perturbations  of  P,Q,R 

Nondimensionalized  perturbations, 
pb/2U0,  qc/2U0,  rb/2Uc 

R  Radius  of  turn 

S  Total  wing  area 

Total  area  of  horizontal  tail 

TAS  True  air  speed 

T  Thrust 

T  Specific  thrust,  T(eng)/pU2d2 , 

where  d  is  propeller  diameter 

U,V,W  Scaler  components  of  Vc 

u,v,w  Perturbation  of  U,V.W 

0,0,0  Nondimensionalized  perturbations, 

u/U  .v/U  ,w/U 
o  o  o 

Vc  Velocity  of  mass  center 

X,Y,Z  Scaler  components  of  aerodynamic 

force  acting  on  mass  center 


lbf  ft 
slug 
lbf  ft 
rad/sec 
rad/sec 


ft/sec 


ft/sec 

ft/sec 


ft/sec 
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x»y»* 

Displacements  along  axes 

ft 

a 

Angle  of  attack  (see  Fig.  1) 

rad 

t 

aw 

Angle  of  attack  of  wing,  a  ♦  iw 

rad 

e 

Sideslip  angle  (see  Fig.  1) 

rad 

6a»*e,6r 

Perturbation  angles  of  ailerons, 
elevator,  and  rudder 

rad 

P 

Air  density 

slug/ft 

8,*»¥ 

Euler  angles 

rad 

Perturbations  of  0,$,V 

rad 

n 

Propeller  RPM 

rad/ sec 

c 

Damping  ratio 

Subscripts 

a 

Aileron 

e 

Elevator 

f 

Vertical  tail 

i 

Inertial 

0 

Denotes  evaluation  at  trim  condition 
(except  for  C,  ,  C.  ) 
o  % 

r 

Rudder 

t 

Horizontal  tail 

T 

Component  due  to  thrust 

w 

Wing 
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1  1 
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(D/l/2pSV^)o 

i 
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C1 
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1 
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INVESTIGATION  OF  STABILITY 
AND  CONTROL  CHARACTERISTICS  OF 
AC 130  LINEAR  MODELS 


I.  Introduction 


Background 

The  weapons  and  fire  control  system  of  the  AC130 
gunship  posses  a  capability  for  extremely  high  accuracy, 
if  the  aircraft  platform  can  be  rigidly  controlled  within 
crucial  tolerances.  A  novel  automatic  control  system, 
termed  the  Sight  Line  Autopilot  (SLAP)  by  its  innovators, 
is  presently  under  development.  The  system  will  be  capable 
of  maintaining  the  aircraft  in  its  firing  orbit  within 
stringent  limits,  allowing  the  pilot  to  devote  his  attention 
to  firing  the  weapons. 

The  replacement  of  the  human  pilot  by  this  control 
system  is  conditional  on  that  system's  capability  to  antici¬ 
pate  accurately  the  aircraft's  perturbed  response  to  a  dis¬ 
turbance  from  trim.  A  comprehensive  mathematical  model  is, 
therefore,  a  necessity.  In  determining  such  a  model,  two 
sets  of  parameters  are  needed  to  describe  adequately  the 
dynamic  idiosyncrasies  of  the  airframe:  stability  derivatives 
and  transfer  functions.  Stability  derivatives,  introduced 
into  the  perturbation  equations,  relate  the  changes  in  the 
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aerodynamic  forces  and  moments  with  changes  in  the  c.ir- 
craft's  state  vectors  when  the  aircraft  is  disturbed  from 
trim.  The  state  vectors  in  this  study  were:  angle  of 
attack,  forward  airspeed,  sideslip  angle,  yaw  angle,  pitch 
angle,  bank  angle,  and  angular  velocities.  Transfer  functions, 
determined  primarily  by  the  stability  derivatives,  are  input/ 
output  ratios  used  within  the  autopilot  to  program  the  re¬ 
quired  changes  back  to  the  trim  condition. 

Purpose  and  Scope 

The  objectives  of  this  study  were  twofold.  First, 
more  comprehensive,  linear  models  of  both  the  AC130A  and 
AC130E  aircraft  were  developed.  The  dynamic  characteristics 
of  the  proposed  AC130E  model  were  then  compared  with  those 
of  an  existing  model.  Second,  data  were  generated  for  the 
AC130A  throughout  an  extensive  flight  envelope  to  establish 
general  trends  and  probable  values  of  stability  derivatives 
and  selected  stability  and  control  parameters.  This  data 
will  be  used  to  compare  model  response  with  flight  test 
results  on  the  AC130A  aircraft.  Aeroelastic  effects  were 
not  included  in  the  model  as  no  test  data  were  available 
and  no  reliable  technique  has  been  developed  for  estimating 
these  phenomena.  Gust  inputs  and  quantitative  determination 
of  transfer  functions  were  not  treated  in  the  study  as  the 
author  feels  that,  until  a  valid  model  has  been  confirmed, 
such  efforts  would  not  produce  useful  results.  Because 
of  the  large  number  of  calculations  involved,  the  iterative 


GAM/AE/72-5 


procedures»  the  dynamic  response  modeling,  and  the  plotting 
of  data  were  all  programmed  in  Extended  Fortran  language  for 

t 

the  CDC  6600  digital  computer  and  the  Calcomp  package. 
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II •  Governing  Equations 

i 

The  equations  used  in  the  mathematical  model  were 
derived  from  Etkin  (Ref  6),  and  include  the  gyroscopic 
coupling  terms.  Assuming  the  airframe  to  be  a  rigid  body 
and  the  earth  axes  to  be  fixed  in  space,  the  complete  set 


of  equations  is: 

Equations  of  Motion 

X+XT-mgsin0  =  m(U+QW-RV)  (1) 

Z4ZT+mgco80cos$  =  m(W+PV-QU)  (2) 

s  QIy+PR(Ix-Iz)+(QR+P)Ixy+(PQ-R)Iyz 

+  <P2-R2)Ixz+hxR"hzP  <3> 

Y+YT+mgcos©sin$  =  m(V+RU-PW)  (4) 

i 

<+«T  *  PIx+QR«z-Iy)-(PQ-R)IXi*(PR-Q>Ixy 

+  (R2-Q2)Iyz-hyR+hzQ  (S> 

«.  *\  «  RIz+PQ<Iy-Ix)  +  (QR-P)Ixz*(PR-Q)Iyz 

♦  «32-pJ>1xy->'xC!+V 


(6) 


GAM/AE/72-5 


Euler  Relationships 


P  «  9  -  i  sin0  (7) 

#  e 

Q  *  0  cost  ♦  V  cosOsint  (8) 

e  i 

R  *  -  0  sint  ♦  V  cos©  cost  (9) 

e 

©  *  Qcost  -  Rsint  (10) 

t  *  (Qsint  ♦  Rcost)tan0  ♦  P  (11) 

e 

V  *  (Qsint  +  Rcost)sec0  (12) 


,  Inertial  Velocity 

z^  a  -Usin©  +  Vcos©  sint  +  Wcos©  cost  (1J 

To  simplify  these  equations  further,  the  following 
assumptions  and  restrictions  were  made: 

1.  The  only  significant  rotors  contributing  to 

gyroscopic  oupling  were  the  propellers.  ■+  h„  *  h„  =  0, 

y  z 

hx  =  ftl. 

2.  The  xz  plane  was  a  plane  of  symmetry.  I  =1  a  o, 

xy  yz 

Vtp  ®  ®  n»T  *  o . 

3.  The  thrust  force  was  parallel  to  the  longitudinal 
axis.  ZT  a  0. 

4.  The  angles  of  attack  and  sideslip  were  not  allowed 
to  exceed  ten  degrees.  The  following  linear  approximations 
were  then  used:  U  »  V  .  a  =  W/U,  3  =  V/U,  W  =  Ua  +  Ua, 

V 

v  =  u$  +  lie. 
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Employing  these  assumptions,  the  equations  became: 


X+J^-mgsine  =  m(U+UQa-UR8)  (la) 

Z+mgcosQcos*  s  m(Ua+Ua+UP0-UQ)  (2a) 

*  QIy+PR(Ix-I2)+(P2-R2)IX2+niR  (3a) 

Y+mgcos0sin*  =  m(U$+U8+UR-UPa)  (4a) 

c?  *  PIx+QR(Iz-Iy)-(PQ+R)Ixz  (5a) 

«  RIz+PQ(Iy-Ix)+(QR-P)Ixz-niQ  (6a) 

Equations  (7)  through  (12)  were  unchanged. 

=  U(-sin0+8cos0sin$+acos0cos$) .  (13a) 


Solving  the  equations  for  time-dependent  parameters, 
the  following  set  of  equations  referred  to  the  trim  conditions 
were  obtained: 

2»Ue  =  P0SU*CX  +  2T0  -  2mg8ineo  ♦  2mU0<R080  -  Q^)  (14) 

2mUo0o  »  P0SU|(CS  -  o0Cx)  -  2a0To  ♦  2mg<oos90oo8*0+a0sineo) 

+  2«UoQo<l*«o)  *  IWVVo1  (15> 

ny«o  =  <>oSUo5Cm  *  JW  W  *  2WRo-Po>  *  2Iy«IRo 

(16) 


60  =  Q0cos^0  -  R0sin$0 


(17) 
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2»U0®o  e  poS^o^y'*®o^*^  *  28oTo  *  *»««>°*Vi,V>o,*"eo) 


•  -  2mU  RU+82)  ♦  2mU  a  (P  +Q  8  )  (18) 

o  o  o  o  o  o  o  o 

2<IVI  -I L)P„  =  PrtSU*b(I  <^+1  C  )  +  2QR<IWI  -!*♦!?_) 
x  z  xz  o  oo  zixzn  o  o  y  z  z  xz 

♦  W«(vvv  ♦  2Vx*81  <19) 

2(1xV4>*0  =  +  2Vo(Ikv4> 

+  Wxz'VVV  +  ^c1^1  (20> 

s  p0  +  (Qosin0o+Rocos^o)tan6o  (21) 

t 

t|>0  =  (Qo8in4»o+Rocos^o)sec0o  (22) 

• 

z.  •  U*(-8in0  +3„cos0  8in$  +o^co80^co8^)  (23) 

1  o  OO  O  OO  0  o 


Thrust  appears  explicitly  in  equations  (14),  (15),  and  (18), 
while  in  equation  (16)  this  term  is  absorbed  into  the  coef¬ 
ficient  Cm*  This  arrangement  was  to  facilitate  an  iteration 
process  used  to  obtain  the  initial  (trim)  conditions. 

Perturbations  about  the  trim  condition  were  then 
introduced.  These  disturbances  were  considered  small  compared 
to  the  steady-state  values  of  the  state  vectors.  The  sines 
and  cosines  of  the  perturbation  angles  were  approximated 
by  the  angles  themselves  and  one  (1),  respectively.  The 
products  and  squares  of  the  perturbation  quantities  were 
considered  negligible  in  comparison  to  the  perturbations 
themselves.  The  following  notation  is  used  in  this  report: 

U  =  U0  +  u,$s$0  +  4>»  etc,  where  the  subscript  denotes  the 
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state  variable  at  its  trim  condition  and  the  lower  case 
character  denotes  a  perturbation.  Assuming  that  inertia 

i  . 

terms  (Ix,  etc)  and  atmospheric  variables  (density,  speed 
of  sound,  etc)  remain  essentially  constant  for  the  small 
time  span  under  consideration,  the  following  perturbation 
equations  were  obtained  (direct  thrust  effects  here  were 
included  in  the  aerodynamic  forces) : 


AX/m  a  u  +  ^a0Q0+R0^0)u  +  +  Urtartq  +  g9cos8 


o'o  o  o' 


-UR$-U8r 
o  o  oo 


(24) 


AZ/m  ®  aQu  +  U0o  +  («0-Q0)u  +  UQa  -  UQq  +  g88in8ocos$0 


+  PoUo®  +  uo®0P  +  6^cos0o8*n<l,o 


(25) 


rtl*  V  *  (Vo-W*Vo>P  +  tIXPo'12Po'2IXZSo*ni)r 


(26) 


A0  a  qco89^  -  rsin$  -  cos6 

O  0  0  0 


(27) 


AY/m  =  6Qu  +  Uo&  +  <R0-P0a0+®0>u  -  Pcu0a  +  g9sin6osin<)0 
+  Uq6  -  OqUqP  +  UQr  -  g4>co80ocos4'o 


(28) 


Wt  •  IXP  -  Ixzr  ♦  <I,VVo-Vo)n(I»-I,lt-IxiV(!S) 
**  ■  V  -  h*  *  ayPo-IxVIX2VM)« 


+  Q  (I  -I  )p  ♦  I  Or 
yo  y  x/p  xzyo 


(30) 


A4>  s  qtan8  sin<|>  ♦  ($  tan©  -P  tan8  +Q  sin4>  +R  cos$  )0 

O  O  O  O  O  O  O  0  0  o 


+  p  +  (tan0o  cos$0)r  +(6  tanO 


(31) 
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=  (qsin^o+0Qotan9osin^o+0Rotan0oco8^o+rco8^o)sec6o 

♦  ♦9os?c0o  (32) 

Equations  (24)  through  (32)  were  then  expanded  by 
linearizing  the  aerodynamic  forces  and  moments  and  expressing 
them  in  terms  of  the  non-dimensional  stability  derivatives. 
The  following  assumptions  were  applied  to  the  linearization 
procedure: 

1.  Cross  derivatives  (Cv  ,  C  ,  etc)  were  considered 

xy  nu 

negligible  (Ref  6:123-4). 

2.  The  thrust  coefficient  was  independent  of  angle  of 
attack  (Ref  6:147). 

3.  The  change  in  side  force  due  to  aileron  deflection 

(C  )  was  negligible  (Ref  3:88). 

0& 

4.  The  airflow  around  the  aircraft  was  assumed  to  be 
quasi-steady.  Unsteady  flow  effects  were  considered  negli¬ 
gible  and  all  derivatives  with  respect  to  accelerations 

•  e 

were  omitted  with  the  exception  of  those  due  to  a.  The  a 
derivatives  were  retained  to  account  for  the  effects  of 
downwash  on  the  horizontal  tail  (Ref  9:4-52;  Ref  4:370). 

The  linearized  aerodynamic  forces  and  moments,  adapted  from 
Ref  6,  w»re  used  in  the  following  form: 

AX  .  |p0SU0(2C  ♦<!„  >u  ♦  KK<Cx  «  ♦  jg-  Cx  4  ♦  C  «.) 

u  o  o  q  fie 

«  ■  ksuo<2Vc*  >“  ♦  ksuo(Cz  • 4 

u  a  o  a  o  q  fie 
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to. -  |ooSU0o<  2 VCBu)u4>0Su|e(C^«^^« 

4Y  *  50osuo<2Cy>ut7posuo<cy6Btrofy,,p*TO^:yj.r+Cy{r5l'> 

^  ■  ^osVaCl)utKsu^(cl88+^lpp+^lrrtclSa<a*cl8r{r) 
All  •  jP0SV(JCn,u+5poSU^<Cn68^n/+1^CnrrtC}{a*atCn{r4r> 


fl 
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III.  Development  of  Model 


Equations 

The  linearized  forces  and  moments  were  substituted 
into  the  perturbation  equations  (24)  through  (32)  and  the 
resulting  set  of  equations  solved  for  the  time  derivatives 
of  the  state  variables.  The  steady- state  (trim)  condition 
of  level ,  turning  flight  was  simulated  by  requiring  the 
following  restrictions: 

1.  All  state  vectors  were  constants  except 
which  was  cyclic,  and  8Q,  which  was  zero. 

2.  Accelerations  of  state  variables  were  zero. 

Applying  these  conditions  to  the  perturbation  equations,  the 
following  set,  referred  to  body  axes,  resulted: 


•  pftsu  p  su2 

u  s  [~S — £(2C  +C  )-ct  Q  ]u  ♦  t-x — 2c  -U  Q  3  a 
2m x  xu  oyo  2m  x  oHoJ 


P^SU  c  .  P„SU  5 

♦  [JSLJELc  ]o  ♦  [A  -2~C  -a  U  ]q 

4m  ^xAJ  4m  x„  o  oJ4 

a  q 

p  SU2 

-  Cgcoseo)8  +  CU0Ro3B  ♦  C-fjr^x  36e 

1  6e 


(33) 
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tl-Trcz434  *  [TEr<JC*+c*u>*of]a  -  Eo;3u  +  c_ferSc^Ja 


♦  Cl  ♦ 


o0sS 

4m  cz. 


3q  -  C^in0ocos$o30 


PASUft 
0  0/. 


-  [Po30  -  E^-coseo8in4>03$  +  — Cz  36e  (34) 

*  0  3 


p„suSJ 


Iyi  ■  ♦  C^fc^Ja  ♦  C^S_cB  3» 


♦  [ 


P„SU  S2 
o  o 


4  ~m 


«-  3q  ♦  «A-W*Vo>» 


p  SU2c 

♦  +  <I2Po-W2IxzVQI)r  <35) 


0  «  (cos^)q  -  (sin$A)r  -  (V»oco80o)^ 


(36) 


i  s  [“^CvfCt,f?,P§"°3u  +  (Po)a  "  Clf8in0osin«o3  9 
m  J  ©  o 


P.SU 


P«Sb 


PGSb 


+  [-2 — 2c  33  +  [  ?■  — C  3p  ♦  C *,*r,”C  -l3r 

LTm  y^  4m  yp  oJ^  »»m  yr 


PftSUft 

O  Or 


♦  C^CO80O  cos$0H  +  C-TU  “Cy6r36r 

II.  C,  C  PASU2b  Clft  Cnft 

tr!Hr:,p  *  posuob[i^ +  i;]u  +  — r~  ci^  *  if36 


(37) 


XZ  2 


♦  [(!*5ui)P  ♦  'r1* 

2  XZ  2  Z 


p  SU  b'  C-  C  I  -I 

♦  [.°°_.<j£+-£i£)+  q  (i+J^-Ji)3p 

4  AXZ  AZ  Z 
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P  SU  b2  C1  Cn  I  -I  I 

♦  i^-rr-^r11  +  +  -  T^)]r 

4  Axz  iz  ;  0  ixz  xz 

p  SU2b  Ci,  Cn,_  p  SU2b  Cl,  Cn- 

♦  +  __Aa]6a  +  _£_2 +  ** 


xz 


2.  u  ,1 

xz 


]6r 


(38) 


I  I 
z  xz- 


C  C. 
n  A  1- 


-T - T^]r  =  P^SU  bt^  +  -iju  + 

lxz  Jx  0  °  Ixz  Xx 


oSU2b  Cna  Cl, 

-S°_ I  _i  +  6)8 

c  +xz  Xx 


♦  [<^r^  -  1!Ro  t(T-^  +  t^)P0  *  r-;^ 


XZ 


X 


xz 


n 


'  I  -I  .  I 


P  SU  b2  Cn  . 

♦  C-S-ifS-tj-E  ♦  _£>  *  Q  <-j — £  ♦  J£)]p 

XZ  ,  X  XZ  X 


psub2  Cn  C1  I  -I 

♦  C-^S-Cr^  ♦  --£)  ♦  Q0(JU  -  1) ]r 
4  xxz  ix  0  iX 

P  su2b  cn,  Cl-  p>  SU2b  Cn, 

+  -2^2-C-t^  +  -^]«a  ♦  S°i**  + 


'<$r- 


]6r 


xz 


xz 

I 


X 


(39) 


<{»  =  ( tan0osin$  Q)q  +  (^ocosOo-Potan6o)  0 


+  (tanG  cosP  )r  +  p 
o  o  r 


(40) 


.  sin4> 

*  =  (^iFi)q  *  (^tan0J  6+  < 


COS<J>f 

cosO 


-)r 


(41) 


As  car  be  seen,  the  numerical  evaluation  of  the  coefficient 
terms  in  these  equations  requires  determination  of  the  trim 
conditions  and  estimation  of  the  application  stability 
derivatives. 
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Initial  (Trim)  Conditions 

The  trim,  equations ,  (14)  through  (23),  were  used  to 
compute  initial  conditions  through  iteration  on  aQ,  T  ,  C^, 

Cm,  and  Cn<  The  iteration  procedure  was  developed  by  person¬ 
nel  in  the  Department  of  Astronautics  and  Computer  Science, 
U.S.  Air  Force  Academy,  and  was  expanded  slightly  for  use 
in  this  study.  The  basis  of  the  iteration  was  the  require¬ 
ment  for  consiant  state  vectors  at  the  trim  condition.  A 
listing  of  the  computer  program  can  be  found  in  Appendix  B. 

Applying  the  restrictions  listed  on  page  11  for  level, 
steady,  turning  flight,  the  following  equations  were  obtained: 


U 


poSUo 

_2_2c  + 

2m  x 


m 


-  gsin6o  -  UoQoao 


(42) 


p  SU 
o  o 

2m 


a  T 


(Cz-aoCx>  - 


0  °  +  rP-(cos8  cos;{>  +a  sin6  ) 
U  o  To  o  o 


mU 


Q  (1+a*) 
o  0 

(43) 

. 

-<}>  sin6 
o  o 

(44) 

• 

ip  cos8  sinb 
o  o  o 

(45) 

• 

)})  cos 6  cos4> 

O  0  o 

(46) 

tan8  /coso 
o  o 

(47) 
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From  Ref  8:203,  the  following  relationships  were  obtained 


for  the  reference  flight  condition: 

..  vc 
"  gtan$ 


(48) 


J  Vc 
*0  s  T 


(49) 


Using  data  from  Ref  1,  approximate  values  of  aQ  and  TQ 
were  computed  to  fcegin  the  iterative  process. 

Iteration  was  first  performed  on  equations  (42)  and 
(43).  It  was  assumed  that,  for  small  disturbances, 


4“  ;  Is  40  +  It  4T 


40  '  to  4“  *  I?  41 


(50) 


(51) 


The  following  equations  were  then  solved  for  Aa  and  AT  by 


Cramer's  Rule. 


du  du 
i  5a  ST 


where  Au  =  UQ  -  u  computed  and  Ac*  a  a0  -  a  computed.  Final 

f  •  • 

values  for  <*0  and  TjQ  were  obtained  by  setting  UQ  «  aQ  =  0 

and  noting  that  requiring  Aa  *  AT  =  0  will  result  in  u 


computed  =  a  computed  =  0. 
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Values  were  thus  obtained  for  a  ,  0^,  P^,  (L,  and  R.  C„ 

o*  o  o  o’  o  x 

and  Cz  were  determined  by  summing  forces  along  the  respective 
axes,  giving 


Cx  8  C?  +  CLsinoi0  -  CDcosao  (52) 

C  *  -C.cosa  -  Cnsina  (53) 

Z  L  O  U  O 

Trim  conditions  were  then  impressed  upon  equations  (16), 

(18),  (19),  and  (20)  and  these  equations  solved  for  C  ,  C  , 

lu  y 

C^,  and  Cn«  A  second  iterative  procedure  was  performed  using 

equations  (16),  (19),  and  (20).  The  matrix  which  was  used 

to  solve  for  final  values  of  C, ,  C„,  and  C  was 

x  m*  n 


• 

Ap 

*  •  1  « 
3p  3p  3p 

5TT  TC 

1  m  n 

si  Jk  JA 

7CT  SET  51; 

1  m  n 

*  « 

ACX 

Aq 

s 

AO 

m 

Ar 

»  « 

3r  3r  3r 

<FC7  TTCT 

i  m  n 

AC 

n 

where  Ap  =  PQ  -  p  computed,  Aq  =  Qq  -  q  computed,  and 

e  •  « 

Ar  *  Rq  -  r  computed. 

Stability  Derivatives 

Formulas  used  in  estimating  stability  derivatives  are 
presented  in  Table  I. 
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Because  the  model  was  to  be  referred  to  the  body  axes, 
the  following^  transformations  (Ref  10:57-8)  were  then  applied 
to  those  stability  derivatives  which  were  originally  refer¬ 
enced  to  the  stability  axes. 


C  S  c  -  <*  c 
m„  m„  o  m 


C.  ~  C.  -  a  C 
la  l  on 


C  —  C  ^  oi  C 
■L  »  o  m 

aB  a  u 


C  +  a_.C, 
n  o  1 


\  *  \  *  °oCs  . 

<1 


C,  “  «  (C,  +C  )  +  a‘C 
Xp  0  Xr  np  0  nr 


C  =  C  +  a  c 
z.  z„  °o  x„ 


c«  +  a_(C,  -C  )  -  a'C., 
np  0  1p  nr  °  Xr 


C  *  C  —  a  0 
X.  x.  “o^z. 

Gg  a  a 


C.  +  a  (C,  -C  )  -  a  C 
1„  o  1  n  '  on 
r  p  r  p 


C  =  C  +  a  C 
2.  ZA  0  xd 
°B 


C  +  o  (C,  +C  )  +  a  C, 
nv,  o  1_  n_  o  1 
r  r  p  p 


C  =  C  -  a  C 
x6efi  x5e  0  zde 


C .  —  C.  —  a  C 

16aB  16a  0  n6a 


C_  8  C,  +  a_C  C  =  a  +  aC. 

z«eB  z6e  0  x6e  nSaB  n6a  0  ^a 


V  =  cy„  '  ®°cy„  %  ’  CV  aocn, 


pB  'P 


C  =  C  +aC  C  =  C  +aC, 

y„  y_  o  y_  n.  ao  lXvi 

rR  ■'r  p  orR  or  6r 


All  other  derivatives  were  either  unchanged  or  were  initially 
defined  in  the  body  axis  system. 
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Model  Characteristics 

After  initial  conditions  and  stability  derivatives  had 
* 

been  evaluated,  equations  (33)  through  (41)  were  arranged 
in  matrix  form 

X  =  A  X 

•  •  •  •  •  • 
where  X  was  the  9x1  matrix  of  time  derivatives  (u,  a,  q,  6,  0, 

•  •  •  • 

Pi  r,  $,  ty),  X  was  the  12x1  matrix  of  state  vectors  and  control 
terms  (u,  <*»  q,  6  ,  $,  p,  r,  <j»,  ,  6a,  6e,  fir)  /  and  A  was 
the  9x12  coefficient  matrix  determined  by  trim  conditions 
and  the  stability  derivatives. 

Mode  parameters  were  estimated  by  factoring  the  coef¬ 
ficient  matrix  A,  using  an  eigenvalue  program  from  Ref  11. 

The  procedure  yielded  one  zero  root,  two  real  roots,  and 
three  complex  pairs.  Factoring  the  complex  roots  gave  the 
natural  (undamped)  frequencies  and  damping  ratios  for  the 
Fhugoid,  Short  Period,  and  Dutch  Roll  modes.  The  real  roots 
determined  the  time  constants  for  the  Rolling  and  Spiral 
modes . 

The  perturbation  equations  were  then  integrated  using 
a  fourth-order  Runge-Kutta  routine.  The  subroutines  were 
arranged  so  that  the  perturbed  response  of  the  state  vectors 
could  be  estimated  from  a  single  control  deflection  or  any 
combination  of  aileron,  elevator,  and  rudder  deflections. 

The  computer  programs  used  to  determine  the  above  character¬ 
istics  are  listed  in  Appendix  D. 
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C130  A  and  E  Model  Differences 

Only  two.  parameters  differ  significantly  between  the 
two  aircraft)  both  due  to  the  fact  that  different  propellers 
are  used.  The  propeller  moment  of  inertia  (I)  and  the  spe¬ 
cific  thrust  <T„)  for  each  aircraft  are  listed  in  Appendix  A. 

C 
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IV.  Results 

Comparison  of  AC130E  Models 

Response  to  a  one- second  pulse  of  0.1  radian  was  in¬ 
vestigated  for  each  of  the  primary  control  surfaces.  Per¬ 
turbations  produced  by  the  model  developed  in  this  study  were 
compared  to  those  produced  by  an  existing  model.  The  trim 
condition  was  turning  flight  at  a  constant  altitude  of 
10,500  feet,  at  28  degrees  of  bank,  at  a  gross  weight  of 
110,000  lb,  and  with  C.G.  at  25%  MAC.  Model  parameters  are 
listed  in  Tables  II  and  III,  and  results  of  the  computer 
simulations  are  presented  in  Appendix  C. 

A  modified  sensitivity  analysis  was  performed  by  equating 
pairs  of  parameters  which  differed  considerably  and  noting 
any  subsequent  change  in  response.  Table  IV  contains  the 
results  of  this  analysis. 

Compared  to  the  existing  model,  the  proposed  model 
showed  an  increase  in  damping  of  phugoid  oscillations.  The  - 
primary  dutch  roll  oscillations  of  the  proposed  model  were 
generally  weaker,  but  the  8  and  r  response  to  aileron 
deflection  showed  a  marked  decrease  in  dutch  roll  damping. 
Pronounced  dutch  roll  was  evident  in  several  of  the  cross 
coupled  perturbations  of  the  proposed  model,  but  the  magni¬ 
tudes  of  these  oscillations  were  quite  small.  The  divergent 
mode  tended  to  be  weaker  in  the  proposed  model.  The  short 
period  oscillations  of  both  models  were  identical. 
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Table  III 

Mode  Parameters  -  AC130E  Models 


Parameter  Present  Model 


Phugoid  Mode 

w  (rad/ sec)  0.146 

n 

0.001 

Short  Period  Mode 

w  (rad/sec)  1.984 

n 

0.636 

Dutch  Roll  Mode 

w_  (rad/sec)  0.866 

n 

0.21S 

Rolling  Mode 

Time  Constant  (sec)  0.663 

Spiral  Mode 

Time  Constant  (sec)  69.31 


Proposed  Model 


0.144 
0.043  , 


2.055 

0.606 


0.968 

0.253 


0.686 


66.09 
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Table  IV.  Perturbations  of  Proposed  AC130E  Model 


Response 

Comparison  with  Existing  Model 

Reason 

a 

Ta 

Less 

dutch  roll. damping 

I 

6 

Ta 

Less 

dutch  roll  damping 

C 

np 

•u  ,  a  ,  q  ,  0 
TalaTaTa 

More  phugoid  damping 

c  ,c  ,c 

V  mu  np  . 

3  >_£.» 

5a  fia  fia  5a 

Less 

dutch  roll  damping 

c 

.  np 

u  ,  e 

Te  Je 

More 

phugoid  damping 

c  ,c 

V  np 

Je  Je  Je 

Less 

l 

dutch  rpll  damping 

'  I  1 

Ve'te 

More 

phugoid  damping  , 

C  jC,  ,c  ,c 

xu  V  nr  np 

More 

divergence 

C1  !  ' 
r 

Tv'Tr 

Less 

dutch  roll1  damping 

I 

u  ,  a  ,  q  ,  6 

5r  6r  5r  6r 

More 

phugoid  damping 

C  ,C 
mu  np 

5r  fir  or  6iv 

More 

dutch  roll  damping 

cn 

Jl>jL 

6r  6r 

l 

P 

Jl>Jl 

6r  6r 

Less 

divergence 

C  ,C,  ,C 
m  1  ’  n 

u  r  p 
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Mode  Parameters  and  Stability  Derivatives,  AC130A 

Selected,  mode  parameters  and  stability  derivatives  were 
calculated  for  the  AC130A  aircraft  in  a  level,  left-hand 
turn  at  30  degrees  of  bank,  altitudes  from  6000  feet  to 
15,000  feet  in  increments  of  1000  feet,  gross  weights  of 
100,000  lb,  110,000  lb,  and  120,000  lb,  and  C.G.  locations  of 
15%  MAC,  25%  MAC  (nominal),  and  30%  MAC.  Results  arc  presented 
graphically 'in  Appendix  D,  From  these  data,  probable  ranges 
of  the  selected  quantities  were  estimated  to  be 


Phugoid  Mode 

w  (rad/sec) - 

n 

£ . 

Short  Period  Mode 

w  (rad/scc) - 

n 

£ . 

Dutch  Roll  Mode 

w  (rad/sec)-; - 

n 

£ . . . 

Rolling  Mode 
time  constant  (sec) 
Spiral  Mode 
time  constant  (sec) 


C 


x 


a 


■  —  0.115 

to 

0.145 

— 0 .025 

to 

0.050 

— 2 . 200 

to 

2.700 

— 0 . 550 

to 

0.650 

— 0 .850 

to 

1.250 

—  0 . 180 

to 

0.300 

— 0.400 

to 

0.900 

■  —  50.00 

to 

70,00 

-0.260 

to 

-0.190 

-0.150 

to 

0.150 

0.040 

to 

0.040 

0.350 

to 

0.7  70 
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» 


-7.700 

to 

-7.500 

-1.780 

to 

-1.730 

-0.040 

to 

0.100 

-3.100 

to 

-2.850 

-9.900 

to 

-9.200 

-0.080 

to 

0.250 

-7.700 

to 

-7.400 

25.000 

to 

-23.500 

-0.090 

to 

-0.060 

-0.564 

to 

-0.540 

-1.644 

(constant) 

-0.749 

to 

-0.740 

-0.050 

to 

-0.045 

0.088 

to 

0.093 

-0.134 

(constant) 

-0.550 

to 

-0.535 

-0.170 

to 

-0.110 

0.388 

to 

0.396 

0.396 

to 

0.420 

-0.173 

to 

-0.155 

-0.075 

to 

—0 • 065 

0.00328 

to 

0.00344 

0.264  (constant) 
0.0264  (constant) 
-0.0920  to  -0.0910 
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V.  Conclusions  and  Recommendations 

i 

The  only  dramatic  differences  between  the  two  AC130E 
models  occurred  in  cross-coupled  responses  which  were  at 
least  one  order  of  magnitude  less  than  their  respective 
primary  perturbations.  It  is  concluded  that  the  existing 
model  is  reasonably  accurate,  based  on  available  data;  how¬ 
ever,  the  qualitative  comparison  provided  the  following 
suggestions: 

1.  Significant  cross-coupling  occurs  in  the 

and  ^  responses,  but  these  are  unaffected  by  the  gyro¬ 
scopic  coupling  due  to  propellers.  On  this  basis,  the 
gyroscopic  coupling  is  negligible;  however,  if  other  cross¬ 
coupling  is  to  be  considered,  these  effects  should  be 
included,  as  the  phenomenon  substantially  increases  the 
dutch  roll  in  several  of  the  smaller  cross-coupled  reactions. 

2.  The  u-derivatives  should  be  more  rigidly  defined 

and  the  contribution  to  Cv>  and  C  due  to  da/du  in  the  tody 

Na  zu 

axis  system,  as  well  as  the  quantitative  change  in  due 
to  axis  transformation  from  stability  to  body  axes,  should  be 
taken  into  account.  These  effects  were  not  included  in  the 
existing  model. 

3.  Flight  test  verification  is  needed  to  determine 

values  for  the  rotary  derivatives.  None  of  the  values  used 

in  either  model  (except  )  can  be  analytically  justified 

P 

as  all  were  obtained  from  empirical  formulae. 
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Further  study  is  necessary  to  provide  information  in 
the  following  areas: 

1.  Further  computer  simulation  using  different 
combinations  of  control  deflections. 

2.  Possible  changes  in  stability  parameters  to  account 
for  aeroelastic  phenomena  and  the  effects  thereof  on  model 
response . 
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Appendix  A 
Basic  Data 

(from  Ref  1,  2,  and  3) 


Horizontal  Vertical 


Wing 

Tail 

Tail 

Units 

Area 

1745. SO 

545.00 

300.00 

ft2 

Span 

132.60 

52.70 

23.10 

ft 

MAC 

13.71 

11.18 

14.82 

ft 

AR 

10.09 

5.02 

1.78 

TR 

C  .52 

0.37 

0.30 

Dihedral 

1.50 

0.00 

N/A 

degrees 

Incidence 

1.50 

-1.75 

N/A 

degrees 

Twist 

-3.00 

0.00 

N/A 

degrees 

Fuselage 


Max  Length 

97.74 

ft 

Frontal  Area 

180.00 

ft2 

Horizontal  Tail  Length 

47 . 33-C ,G. (13 . 71) 

ft 

Vertical  Tail  Length 

46. 03-C.G. (13.71) 

ft 

Propeller 

RPM 

106.80 

rad/sec 

Diameter  A  Model 

16.00 

ft 

E  Model 

13.50 

ft 

I  A  Model 

139.00 

slug/ft? 

E  Model 

175.40 

slug/ft^ 

Tc  *  0. 375-0. 01A-(Q. 00338-0. 0009A)AV 

(a) 

T  *  0. 275-0. 008(A-10>- 

-CO .  00246 

l-0.00007(A-10)3AV 

0-) 

where  (a)  is  valid  for  A  <  10000  ft,  (b)  is  valid  for 

10000  ft  <  A  <  20000  ft,  and  V  =  TAS-120  (knots). 
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Z  axis 


Y  axis 


Figure  1,  Eody  Axis  System 
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Appendix  B 


C130E  Model  Response  to  Control  Deflection 
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Appendix  C 


Stability  Parameters  for  the  AC130A  Aircraft 


l 


i 


I 
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VARIATION  WITH  UCIOHT  AMO  ALTITUOC 


VARIATION  KITH  C.O.t 


VARIATION  WITH  UCIOHT  AND  ALTITUDE 


Vn 


VARIATION  WITH  C.O.i 


Picure  9 ,  Short  Period  Parameter 
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VARIATION  WITH  WEI DMT  ANO  ALTITUDE 


Iff. 00  170,00 


VARIATION  WITH  C.O.  POSITION*  NONE 


VARIATION  WITH  WEIGHT  AND  ALTITUDE 


«• 


VARIATION  WITH  C.O.  POSITION*  NONE 


TAU(SEC)  TAUtSEC) 

t.4f  Ml  Ml  t.W  Jf.n  W.W  •••TO  10.TO  iO.m 


Spiral  I.:ode 

VARIATION  KITH  WEIGHT  AND  ALTITUDE 


VARIATION  WITH  C.C.  POSITION!  NONE 


Rolling  Kode 

VARIATION  WITH  WEIOHT  AND  ALTITUOE 


VARIATION  WITH  C.C.  POSITION!  NONE 


Figure  11,  Spiral  and  Rolling  lodes  Parameters 


VARIATION  M(Tn  UCIDMT  UNO  ALTITUDE 
(NO  VARIATION  HUM  C.O.  I 


Figure  14.  Stability  Derivative  Variation  (Cx 


Figure  15.  Stability  Derivative  Variation  (C«  ,Cg  ) 

"4  Q 


Figure  16.  Stability  Derivative  Variation  (C^  ,Cm.) 


vmihtioh  kith  kcibht  on  hltituoc 


v 


VMtMTXO*  KITH  C.0.1 


CH4*  «  -I. Ml  HMD  UM  CIMTHHT  nt  HU  HLIOHT  CMOiTIOKI 


Figure  17.  'I  ability  Derivative  Variation  (0y,  ,C», ,  V. ) 

%  na» 


53 


Figure  13.  Stability  Derivative  Variation  (Cy^,C^ *Cn^) 


w 


j  VARIATION  U  E  TM  MEtOMT  AMO  ALTITUDE 


VARIATION  WITH  UEIOHT  AND  ALTITUDE 
INO  VARIATION  WITH  C.O. I 


VARIATION  KITH  WEIONT  AND  ALTITUDE 


Ficurc  20.  Stability  Lorivative  Variation  (Clc,Cn.  ,Cm 


*  /  ^*1  |0m  |C 

Figure  21.  Stability  Derivative  Variation  p  p 

\  Cy^  >r'l& 
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Appendix  D 


Computer  Programs 


p* c**'AM  c,,,-'PUTrs  r;iT:u  rcrniTieNS#MCCE  rapametfpn 

AND  P£CT,P'»:*'  =‘Z*;F')N'i£tP’Ir.'r  ST AT£rtiNTc  NOT  muCEC 

ESTTMATf:  -?:iCcv  ANO  4JI  SEOONOS 

THIS  PROGRAM  Tc  SET  UP  Prt°  TM*  c  NO^EL 

IF  A  MOOtl  I «  UaET*  CHANGE  FOLLOWING  terms... 

Tr=T/<*i*i.:;?> ,  GT=i.c‘f2*T<',  1=120.: 

DIMENSION  Xi.  ( 2^15)  ,X?(?S«>  ,  v 3  (  ?5r») ,  X4  (2E5 )  ,X5  (25S)  , 

*V6  C255)  .X?  (?.>«)  ,XS  <2fr)  ,X9(25  5>  ,X1  )  (255)  . 

*Xii('*5E)  |X12(?55) 
ni^-NSICN  t  i  {  P5«=)  ,y(!2) 

CCMMCN  AA(c,i2) ,v(25?f 12) 

DIMENSION  KM9) » 7FTA (°) »P0<9) »THALF(9)>A°EAL(2) 

°EflL  K,K1»  KEtK3, Kfc,KP  tKF*  «7  ,K8 ,  «9.  K 1  ,  Ki  1  ,K12  ,K13 ,  K14 , 
*K13,Klfi,*l7f<l«f<2D»K21f''22»K2?,K2.<f  K?k»<26,K27,K28,K29 
PEAL  I,!XX,IYY,nZ,:XZ,M,Ml,LT,LF 
INPUT  CftPaMCiros  UHIO”  ANr  CONSTANT 
READ  2t  J,S»c»?,C  iAo,ei,Rr'K,AP 
INPUT  PARAMETERS  ^ICH  VARY  WITH  FLIGHT  CCNDITTCN 
RPAQ  l-1fUfALTC»cMlffLAcS|W.TPMP 
WHERE  U  IS  TAS  IN  *NOT*fALTO  IS  oR^SS  ALTITUDE 
IN  FEET»phI  IS  3f N<  ANGLE  1M  OtG?EFS,W  IS  GROSS 
HEIGHT  in  L°,ANC  HP*  IS  Fr«-r  STHA*  AI°  TEMP 
READ  lr5,IXX,IYY,I7Z,iyZtItCG 
COMPUTE  PA=AMET'E$  NEE CFO  TO  DETERMINE  AERCCYNAMIC  DATA 
°PM=1 32  3 .D*PI/3C.P 

PHI  =  omj/57, 2957fl 
m=W/o 

*HOa.  0.2376**  <1.0-.,.  I  GO  068825*  ALTO)  **4.255 
f=plafs/ioc.? 

DEIV#U-120.15 
UsU^jf},  0/2603.'’ 

VSOUNQs49.1*SORT(TEMP) 

Ml  =  tl/V$OUND 
ALTO-AITO/IojO.O 
CO^sl.l*(.'J?6+,0Ci2*F) 

Ks.«n2-.b2e*F 
CLMIM0G=.2S+.?*F 
COMPUTE  CONSTANTS 

Ki=Ru 0*3 *U **2/2.9 
K2=K1/M 

K3=owo*S*U/(2.0*M) 

K4=K1*P 

INITIAL  VALUES  OF  A  ANC  T  TO  PFGIN  ITE°ATION 
IF  (PHI.pc.C.O)  GO  TO  4A 
Pls’J**?/  (G*TAN  (PHI) ) 

PSIOCT  s'J/Rl 
GO  *“0  47 

46  PSin0T=3.0 

47  CONTINUE 

IF  ( ALTO.GT.lO.u)  GO  TO  21 
TC=.375-.01*ALTO-(.ur33P-.u0; 39*ALTO)*DELV 
GO  TO  2? 

21  TCs.27S-.j?P*(ALTC«lr.0> 

*-(. j"24^~. CC306G*  (ALTO-) -.) )*DELV 
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. . . 


2’  CONTINLP 
rT- J.dZ*,*TC 
Ts-cr^vi 

CLLs'i/tKl^CfStrHT)) 

CLL  1£*TC4.9:*F 

IF  <*ULT.W#225>  30  Tn  7 

^LJ-<C*74,*<,254Tr<-2«C>p)*U.f,  +  .?33f‘11-.2>5)) 

3  CLL  As6  •  3  +  ,  4*  20*T  O-t-2.  p  *f 
4  AsOLL-CLU  )VCLLa 

ITERATION  CN  ALPHA  ANC  THRUST  USING  CRAMERS  RULE 
DO  1 J  J  =  1  f5  ” 

IF  (J.LT.2)  00  Tft  14 
CLL^.PS+.lE^r-^g^p 
Ir  (Ml«LT,  C.22R)  GO  TO  1«? 
rhLTr(^*3t*<‘,5*Tr+?*  J*»3?2(Ml-.225> ) 

oU  Tr  34 

15  CLLA  =  S»  7+,  4  2:>*TC+2«  O*!1' 

14  cll=cll:+clla*a 

C0*C00  ■*•<* (CLL-CLMINDO  **2f  ,  00  5 
f  TsT/Kl 

00*03*  <1.+*1a*F*CT) 

OOA*?, *<*CLLA* (CLL-CLMINOR) 

CX* CLL* SIN  (A) -CG*COS (A) 

CXA= (CLL-CCA) *OCS ( A) ♦ (CD+CLLA ) *SIN (A ) 
C7=-rLL*C0S  (A)  «CE*SIN(A> 

CZA,=  (CLL-CCA)  *STN  ( A ) - (CO+CLLA > *COS ( A > 

THaATAN(TAK(A)«nC3(°HI)t 

cs-?SICCT*SIV'(TH) 

G=°STCor*cC? (TH)* STM (PHI) 

PsPSTDOT^CCStTH) *rcs (°MT) 

MOOTs X?* Cx  +  T/f1-v>*SIN  (TH)  »U*0*A 
AOOTsK?#  (C2“A*CX)-(A*T)/cm*(J)  ( 1 ,  +  A*  *2) 

*+  (G/*J)  *  (COS  (TH)  *CCS(C«I  )*A*$IN(TH)  ) 

CTHOAsCOS (PHI) /(l. > (SIN(3HT) *SIN( A) )**9) 
opjAr-rsncT*co5- (th)  *othoa 
OQDAs-PSnCT*STN  (rHT ) *SI” (TH)  *OTWOA 
OROAs.psxncTvccS(F:«I)  *ST^(TH)  *OTHOA 
^°^!j,,!^C’r’4CnS(TH>¥0TH':J4-u*(Qf  A*00PA) 

^f?^r(CM'Cr"A,,CXfi,*-T/(M^)  ♦  (G/U)*  (STN(TH) 

^OOO^A^r^”  1A*’SIN(TH,¥C0S<RHIMCTH0A)+?*'?4A*( 

OADDTs-A/ 

OCTsOUOOA*CA33T"OunOT»OAGOA 
OEL  As  (Cl,’,inMOOT-riAr)n*lir>QT  )/0ET 

0ELT*(CAU0A*U30T-rUDrA*Ai)0T)7CET 

AsAtOELA 

TsTfOELT 

TC=T/(K1*.8?6)  ' 

IP  CONTINUE 

TRANSFORM  INERTIA  TE°MS  TO  *>QOY  AXIS 
ALPHA2A"{37.?0'>7A 
IXZsTXZ*(l,-ALrHA/4»7‘3) 
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IXXSslXX 

IZZS=IZ7 

ixzssixz 

IXXsIZZ3*SIN<A) **2*iyx?*C0?(A)**?+?.*IV7?*SIN{A) ♦roSfA) 
IZ?=IXy3*STMA>  ♦*?fII7S*C0r:(A)**2-2«*IX7‘r*STMA)*''C<?fA) 
^  «  ZxZsI*zSMCCS<A>**?-SIN'<A)**2>MIZ7S-IXXb>*SINtA>*C9S<A> 
COMPUTE1  CONSTANT? 

K5*IXZ/IZZ 

K6*IYY-IZZ+(TX7**2/IZZ> 

K7sixz*(i#:+(ixx-iyy)/tz?) 

K8sK5*I*9pv*t*.C 

K9*K1»C?AR/IYY 

KiOs(I?Z-TXX)/TYY 

KUalXZ/TYY 

Ki2*-4.-*I*PPM/IYY 

KIZsIXZ/TXX 

K14aIX7*((lYY-IZZ)/lXV-l.o 

Ki5=IXX-TYY+(lXZ»*2/IXY) 

K16*T*PPM*  it,. 

K17*TXX-(IXZ**2/!Z7> 

XlfisIZZ- (IX2**2/IXX) 

K2i.sK5*Ki44K5 
K21sY5,Klt>*^7 
K22=K5*K154IO 
K23=X4*  (1.  :-*5*X!?> 

K?4s<6^K1‘?4Y14 
K25aK7*<i3+«15 
K26=K3*<13**16 
K27  =  ’<  11/X9 

K28a-Kl^/XC 
K29a-Ki2/K5 
TC*T/<K1*. 575) 

CLL0s.25  +  ,U*T''  +  .a;*P 

IP  <M1.LT,C.725)  r,o  t0  ^ 

CLLAs(6.7>,475*Tr*2,i*r)*(i,^4,733(^i-.??5)) 

60  TO  17 

16  CllA=6.?4.4iT>*TC  +  2.?*F 

17  CLLsCLL^CLLA*® 

CO*C03  +  KMCLL-ri»'IN0C)**24t  j"5 
CTaT/Kl 

C0sC3*(1.4.13*P*CT) 

COA  =  ?.*K*CUA* <CLL-ClMIUOG) 

THsATAM(TAMA)  *CC«(f‘HI) ) 

Ps-oci:.OT*Sri(Tri) 

0=PS1H0T4CCC(TH)  * «?  !►!  f  PH  I ) 

R*P3IC0T*CC5(TH)  *rn<;(r>Hi) 

PYs(C.o»A.  (G/U)*,'C^(Tr4)  ♦«?I*|(r  HI)  )  /X7 

CMsX?7* 

rNs(X2‘.  +  *.*.)  /K?i 

ITC  *  AT  ICM  ON  ^  ,  C I  *  C “  031*^  r  *■  A M c. ? ~  PMLc 
nn  4'*  jj=i,r 

r  D  ^  T  ^  L  ♦<  3  *  *“  “ )«•>•',  +  '<  7  *  =  *  4  <  5  *  C. )  /  "*  1  7 

COOT  =  Kc¥Cl  +  Yiv*u,t^4XJi#(>‘,*,f2“P<1*?)4X12*5! 
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oH0*s  («'4*<0*Ki  i*rjL I  ♦  KlA*G*R*'<i«*P*C4f''lS'M,'>  /IMa 

''P  )*><**»-  -  b  ~  3  " 

0PC'>,'L=M4/K17 
CP33PMr<,*K5/Klft 
C0DTCMsK9 
,  roDccur;.j 
COOO^Ns?  ,  <3 
DRDOCM=  *■ .  A 
CPinCLsK^^KU/Kl? 

0RD0CV=K4/^lH 

OFT  lsQCQfCt'*(QPO'?CL*nPDCf'N-DR3DCL*CcCrCN) 

DFLCv*C?OT*  (JPrr.CN*0?r^r*ri  -DPOOCL’CRFCr^)  /pptj 

nFtC»l={)OnnrM*  ( PC0T*DrDP^f#-R00T*0oDDCL>/!5rTl 

D£LCL=CCC^rP*(KCCT*3f!n:CL-FD0T*npb3CN)/1cTi 

CM=C“+nPLCw 

.CUCUCCLCL 

CNsCN  +  nrLCK* 

«  CONTINUE  -''-SSSKsa 

MISCELLANEOUS  PARAMETER?  Vw.vjjJSsS-*'’’^ 

awr=alp«a  +  i.p3  VS?-5' 

CALL  ETA  (CLL»  TC>  Ff  FTATf IL) 

CL  A  T  AIL  =  3, 2-"  S*  FT  AT  AIL 
CLAWs(.li4+.,52*F) 

CLW=.2t+l. ^2*F+TLAM*AW« 

COtW=.v.:-52+.:')4A*F 
IF  (TC.GT.j.l)  p*C  TO  ?3 
OCMCLL*-.  24'»  +  .  3?  *TC*  ( .  f 22-.72*TC)  *e 
GO  TP  24 

ll  PONTINUr*2C^*r^<TC”,1>*<*f’5+'ie’(Tr**in#F 

CLAW=CLAW*F7«29R,,8 

nETAlsSQRT(1.0-Ml**?) 

yxis(AP»Ml**i')/(2.3*(A°*P£TAl42»C)*BETAi) 

LT*«*7.3‘»-CG*C°AR 

LFsLT-1.2? 

COMPUTE  DERIVATIVES  REFERENCED  TO  STAQTLTTY  AXCS 
CTUs-3.G*CT 
CXU=CTl’*COS  (A) 

C2U--(CLL*Mi**2)/(i,.M1**2).CTu#STN(fl) 

CMU* JCMT^LL^tCLL^Ml^*?/ 

CXAsCLL-COA 

CZA=-CLLA-CC 

CMAsCCMOCLl*OLLA 

CXQsi.t 

rZO=-?.?*ETATAlL*  (LT/C°AR) 

CHQs-2,2*FTATAIL* <LT/CpAp>«*2 
CXAOsO.a 

CZA0=-2, C^ETATAIL^OEFJA* (LT/C°AR) 
CMA0s«2tj,fETATAlL<lCFR0A<1’(LT/C?AR),|t*2 
CMOE*-. J2a?*&Z.29F7S 
CZrErMC°AC/LT)  *CTATAIL*CMOE 
CXDE=-.?BD 

rYt1s-  ( •  j!55-,  0  36*  Ft  ( ,  33:-,  35.1 8* F)  *  TC)  *57, 29F78 
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CL3-- (•?'  "of-,  r  •;  c  ” 

^Nns(,r. ’1374.;*  ^♦e.,fjlc;»Tr;'j  *57#,3,:7g 
rF  (TC.GE.  ...f)  GC  TO  5  “ 

CY'»^=-.6i.5*...c5*P  +  ( i.«.7i4-2.  5*F)*TC-<1 .?F-2.C3*F)  Mr**2 

CL^Tr-,  0*35-, uC?4*e- (.022*.  i)9*F}*T0*  (• C9*95*F) *TC**9 
•GO  Tn.  6 

*  5  CY9T*-.<«jO-.63«*e 

CL9Tr-,  JAI4  3-,  j11c*?4  (.LC.q-,  026*F)  MTC-.?'5) 

6  C»inTs(.  -j.3  5-,0  jCt*FX.2J?*F-.  JO  75)  *TC)  *5  ?  .  29578 
CY -R  = « . ^ 46*5 7« 2 9574 

‘  CLDRM  •  D0  346-( .  CC3  21/A.  C>  *  ALPHA)  *57,s>o«57A 

CNOP=-. 3„i£*97.29«78 
CYOAsJ. 3 

IP  (CLL.G.T.1.3)  GO  to  11 
CLOAs-. :6j2-. 01?*P+, t29» (CLL-.80) 

GOT^l? 

11  IP  (rLL.GT.1.2; )  PQ  TO  it 

CLOA--,  j544-.  0l72*Ft  (  .3143+ .93?*F)  MCLL- 1.30 
GO  TO  13 

12  CLDA=-,  3K16-.dll4*F+,029MCLL-i«2O 

13  CONTINUE 

IF  (ALcHA.GT.2.u)  GO  TO  1 

CND A= (. w rc 063+ . 0wG02*ALPH A) *5 7.29578 

GO  TO  2 

1  CMDA=(,  J  ’j  )[  57*  .  C’H  012F*  ( ALPHA-2,  c)  )*57. 29978 

2  rv°=-. ulf7*CLAK-°. *A*CYnT*LF/9 
GLPS-.CSJJ 

8  CNPs-. t9*CLK-. C S15-2, *A*CN3T*LF/B 
CYPs2« f  *CN°T 

CLP  =  ..?0*CLW*(l.fXXl)  ♦,  3  51<)*F-CL3-CY9*7*CCS(A  )**2/P 

M2.D*CN9T*CL3T/gy*T  C* 

CNRs**,  ,,2*0 LK**?-,?*CDwW+c.  u*CNBT**2/CY9T 
CORRECTION  FGR  MflCH  EFpFOTS 

IF  (M1.LE.  .225)  rc  TC  47u 

CYUsCY~Ml  •  ♦•  15MMi-.225) ) 

CNL^CNP*  (1  .*.2*  (•'1-.225) > 

CYnT=CY9T*(1.+,2^(Mi..225)) 

CN0T=CK'!?t*  ( 1,  +  ,4*  (Ml-,  2?R) ) 

CLp5CLf--.t2MMl-.225) 

47 9  CONTINUE 

TRANSFORM  PORIVATIVES  TO  900Y  AXES 
CXsCT+f*CLL-CO 
CZa- (CLL+AMD) 

CXA^sCLL-CDA+A*(Cn+CLLA) 

CZA35A*(CLL-CCA)-C0-CLLA 

CMA1sCHA+A*CMU 

CXU?=Cyj4CZf*A**?-A»  (CXA  +  CZU) 
r7U?5C7'J-CXA*A**2  +  A*(CXU-CZA) 

CMUr)sCMJ-A*CMA  • 

CXCnsCX^-A  *C7Q 
CZC  1  =  C7C *A*fXQ 
CHOTsCPO 

CXAjr,5CXAO-A*C7AO 
C?A3‘  =  C7A!HA*CXAr 


Cvn£~=f  xa£-A*c,',r  ' 

CZCi£1=C?3=  +  fi*C)ar 
CHD2!’='“nr' 

CY-‘?=CY° 

CL«^sCt!3-A*C‘if> 

CNnosCNT+A^rtn 

rYP9=CYP-A*CYR 

CLP°=CIP-A*  (£L°+CNC>  +^NP*A**2 
CNP>»=Cf.=-A*  (/'NK-0Lc)  -CLB*A**? 

CYRSsCY^+A^rYP 

CL<?3=CIP-A* (CNn-Cl°)-CNP*A**2 

CNkq  =  CNR4-Ay  (CL«?+^NC)  +clp#a**? 

CYOABsCYOA 

CYns=»=CYiP 

CL04R=ClCA-£*CNDA 

CLD?nsCt3P-A*CNn^ 

CND AcsCNuA  +  A*CL" A 
CNn9o-CM3(H  A*CtCR 
COMPUTE  COCFFICIFKT  matrix 
U0=U 
AC=A 
00?;/ 

THOaTH 

pr=° 

>50  =  '? 

PHI 1=°HI 

AAlsPHC*S*UC/(2.*M> 

AA2=AAl*U(J 
AA3= AAi/UO 

AA4=PHC***CFa»/<4,*M> 

aa«;=>?ho*s*uc*c?A5’/(2,»iyy) 

AA6=AA5*U0 

AA7=AAF*CRAP/?. 

AA8  =  3HG*5*«/{<».*M) 

AA9=RP0*^*U:*e 
AA1 J=AA9*U j/2* 

AAll=AA9*n/4. 

AA12=1.-AAA*vJZA'?0 

AA13=IXX/IX7-IYZ/IZ7 

AA14*IZ7/IXZ-IX7/IXX 

AAlisRHOTS^C’APTHC^CXAnp/t^.VM) 

AAiSsAAi^/A A12 

AA17sl,+AAlF*AC/U0 

AA18  =  RPO*S*li?*r°AR/(ti,*M) 

AA19=i.-AA4VCZA'1B-Ar*CXADs<) 

AA<l,6)sAA(1,7)=AA(l,5)=7A(l,q)=AA(l,10)sA«(l,12)a0,G 
AA  (2,6)  =  AA (?,7)=AA(2,9)  =  AA(2,1Q)sAA (2,12)=a,e 

AA(3,9)  =  AA  (3»1D=AA(c»12)  =  ?»n 
AA(4,1)=AA(u,2) =AA(4,4) =AA(4,5)=AA(4,6)s9,o 
AA(4,P)  =  AA<4fl.'  =  AA(t,H)sAA(4,12)sO.C 
AA(«i,q>  =  aA(5,,3)=AA(5,il)=3,9 
AA(6,2)=AA(f  «4>  =  AA<S4>=AA<6t9>  =  AA(6»H>su.C 
AA(7,2)=AA  (7,4)=AA(7,9)=aA (7,9)=aa (7tli)=C.f 
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AAH,l)  =  AA(c,2)sAAn,r)=AA(<),  •)=0.  3 
a  a  ( •) ,  g )  a  a  a  ( r ,  1  * » = m  c  > , : :  >  =  1 4  ( « ,  1 2 )  a : .  f 

A£(Sfl)slA(c9?|s*.ftrot>:)S;A(iys)Sf>Ar9v(|, 
ftA<  9, <3)  si; (f,  1  »-* A  (<••,!:>  =  4 u  3,i2)  si.f 

AAfl,l|s(lfiI*(9t*cx*Cv,,p)-A;*0‘. 

,  **AAi«:*  (AA**  <‘,,#'*7*f'7t.  v»)  j/U. ) )  /  AA17 

AA  <1  »2)  s  <  AAt*C*A°-'l.'«  *r. A  is*  Ml*C2A  2)  /  A A 17 

A«  <l»3)s  (-AC*U:*AA16*  (Cv^G  +  CX  AD"*  { 1 , +Afc!**rZpa) /AA  12) ) 
*/A  A 1 7 

AA  C 1 1 1» )  r  C-i*C0S(TM3)  -  AA16*(6*SlN(Ti))<,C0''(«5HlD)/U!‘)) 
*/AAl 7 


AA  (1,5)  =  (Uf‘*R3-AA16*C'?)  /AA17 
AA(l,8)=AAi6*5*rOS(TH.)  *$  IN (PM  I**)  / ( U  j»AA  17) 

AA(2,1)  =  (AA?*  (J.^C^+CZ'.^-A1.*  (2,*rXfCXl">) ) 

♦♦Oj*  (l«  +  A"**?)/ti')/AflQ 
AA(2,2)s(AAl*{CZ<>P-Ar*CyAD)+Aj*''0)/A4jC 

AA  (2,?)s  Cl  *4  4?**2*AA4*  (CZ3  !-A'»*CyCa)  )/4A19 
AA(29<i)s((G/,.i?)MA,*''QS  (THw)-SIS(TH0)*CO,>(PHIQ)  ))/AAig 
AA(2,5)s(-P£-AC*5?)/AA1?  ^ 

AA(?,3)a«(  (G/Uw  )  *COS  (THC  )**?I:(  (PHi:)  ) /AA19 
A  A  n,l)sAA5*<2.*CK*CN,,»',)*AA(2,i)*Aa**CMA0»> 

AA(?f2)  =  A36*,JMA°  +  /,4(rf?)*AA7*CMflr'!» 

AA(3,7)  =  AA7*(ruC'-fAA(?,?)*CMADP) 

AA  n,4)rAA  (?,4)*AA7*rMAQn 
AA(’»5)«AA  (?,5)  *AA7*fMAC» 


AA(3,G)s  (  (IZ7-IXX)  ♦■5C-2.*IXZ*P0)/IYY 
A  A  (1,7) a <  ( IZ’-IXX)  *»,.  +  ?,  *IXZ*SC-»?CX*  1*4 . )  /I  YY 

AA(3,3)  =  AA(?,3)*AA',*rvAD0 


AA<4,7)=CCc(«»HI.{*) 


AA(At7)a-SU(PPI'>) 

AA(A,  8)  a-PS^COT^CO*)  (TH^  ) 
AA^jDaPHC^CY/M-^O/Uw  +  A^PS/UO 
AA  (*!  ,2  )  a°3 

aa(5,a)s«g*stn  (TH:)*cpi(?Hn)/ir 
AA(5,5)aA Al^nY^g 


AA  (5,6)sAA9*CYpn  +  A0 
AA  (5 ,7) bAa«*CYR?-1, 

AA (5 , 8 ) s  +  G4rOS (TH  c ) *COS ( PHI  3) /U° 

AA(6,l)sAAc* (CL/lX7trN/I7?)/AA17 
AA  (*>,3)s((  (IXX-TYY)/lZZ  +  l,)*oi 
*+(  (IYY-IZ7)  /TXZ-IXZ/I77)  *"‘‘*gPK*I*A./IZZ)  /AAl  J 
AA(<>,5)  aAAlc  *  <CLr{?/iyZ  +  rN03/IZZ)/AAl? 
AA^SlatAAlt*  (Cl.P<3/7X7  +  r‘*JPP/IZZ) 
♦♦QC*(1**(I>X-IYY)/IZ7))/AA13 
AA(6,7)s(AAll*(ClP°/iyz+CNR3/I72) 

♦♦QvJ*(  (IYY-IZ?) /IYZ-iyZ/IZZ)  >/AAl 3 
AA<7,l)sflAO#(CN/Tx2fCl/IXX) /AA14 
AA(7,?)r ( (<IXX-IYY>/IX7+JX?/IXX)*PQ 
*  +  ((lYY-I2Z)/IXX-l,)*rr'C*t-'rM*I*A*/lX7)/AAlA 
AA(7,5)sAAlC* (CNpn/IXZ+^«  )g/TxX)/AA14 
A A (7,5) a( A  All* (CNPVIXT  +  ClPg/IXX) 
•♦crMiy7/ix»<m-m)/IXZ))/Afll4 
A  A  (7, 7>  a  (  4  A 11*  (GNRg/IYZ  +  rt.pg/IXX) 

*fQ0*< (lYY-IZZ)/IXy-l.))/4Al4 
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flA(<J,‘»)=TftN  <TW  )  *ctm  (pht  '  ) 
AA(Sf4)=F5I?1T/rrS(7Mj) -P.?*TAN(TH:> 

AA(D, 7)=TAN  (TH,)*nrS(Pwy»') 
*AA(9,?)=P!N(°HT:)/C0S<tHl) 

AA<S,,4)  =  PSIC'>T*T5M(THr) 

AA(9,7)sO0P(PHT')/r9f(THn 

AA(l,ll)  =  (AA2*CX''r3'fAA16*AAl*CZr)Ea)  /AM7 

AA(2,  11)  =  (AAi*  (r.,C'~*-AC*PX,'E,'))/aAi9 

AA(3,ii) siA6*Cw9H?+AA7*CMA9a*AA(2f il) 

A  A  (5  , lb )  * C  •  i 
AA(5,12)  rAAj»CV?»n 

AA(e,.m  =  ((CLOAG/lX?+CNDA3/IZX)*AAi3)/AA13 
AACe^ZJsUruDPS/IXZ+CNnce/TZZJ^AAl  j)/AAl? 
AA(7,lt)r(  <CNQA',/IX7  +  CLDA;i/IXX)  *4A10)/AA14 
AA(7,U)  =  <<C'!DPr,/lXZ+CL0‘  3/IXY>*AUj)/AA14 
CALL  Mrr,vAt  (9,AA,WN,7ffTA,P0,TwALF»ActAL) 
WHERr  WN  IS  NATURAL  FRrPU£NCY , Zlt A  IS  DAMPING 
RATIO, PD  IS  PE?IOr,TWALF  IF  TI«E  TO  DAMP  TO  HALF 
AMPLITUDE  ANC  A7EAL  IS  MATRIX  OF  TIME 
CONSTANTS  TCP  POLLING  ANT  SPIRAL  MODES 
INT'EG0 ATION  FY  RUNrE  K'ITT A 
THIS  ARRANGEMENT  RFCCROS  EVERY  10TH  VALUE 
SET  INITIAL  VALUES  C?  cf RTUR°ATI0N3  TO  7£P0 
00  454  KJri,i2 
454  XC1,KJ)=0.C 
TK1XUC 
T=D. 

VALUE  CF  OELT  USE^  WOULD  PE  SPECIFIED  HERE 
11=1 

DO  314  J=l, 12 
314  Y(J)=X(1,J) 

318  PONTINUE 

DO  316  Jsl , l’ 

316  CALL  RUNGE1  <T , MFTOUNT, TELT , Y> 
n=im 

DO  317  J=l,12 

317  X(II,J)=Y(J) 

T1 (II) =T 

IF  (II.Lii.2FJ)  GO  TO  318 
JJJ=II 

CONVERT  PARIAN  MEASURE  TO  CCGREES 
no  99  JJ=1,JJJ 
TII*JJ 

X1(III)=X(JJ,1> 

X2(III)  =  X(Jj,2)*!i7, 29578 
X3<ni)=X(Jj,3)*‘‘7.2957P> 

X4(III)  =  X(Jv!,4)*57. 25578 
XSaiTJsYtjj.BJ^cy^cjTe 
X6(IIT)=X(JJ,6)*:>7. 29578 
X7(III)=X(JJ,7)*57. 29578 
X8(III)=Y(jj,8) *97.20578 
X9(IlI)sY(JJ,9>*57. 29578 
XU  (III)  =X(uJ,  ID)  *57.29578 


Xli(TTI)*XijJ,li)*37,r5r?4 
X12(III)  -5X(vJ,12)  *‘J7. 29*7  4 
ca  CONTINUE 
SU0:?CU*TIN£S 
TAIL  EFFICACY 

SUG9CUTIN*  eta  (rLl»TCiF»ET*TAIL)  vy 
IF  (TC.GT.C.lu)  GO  TO  2 
TF  (CLL.GT.1.2U  r.q  TO  1 

fTATail  =  1.  j  +  fl.S^-.S^F-Mi.’S-.^^F)  *(CU-:.3C)  )*Tr 
GO  to  4 

1  ETATAlLsi.:4(i,7E-.7*F+(.5Gfi.5*F>*(CU-i.2C))*  (TC-u.iO) 
GO  TC  4 

2  IF  (rLL, GT •!.?:>  GC  to  3 
FTATAILsl.l?>-."c:*F4(.i?r-.';G4F)*4CLL-u.«'?) 

X4<.Oi-i.iG*F+(*7754«c*F)*<CLL-C»P3) ) *TC 
GO  T*  L 


2  ETAT*.!Lrl.l7;>-..,7*F4(.i625  +  .i25D*F)*(rLL-1.2(') 

X4(l,  08-.96*F4<.?i294.a75w*F>MCLL-1.2C>>  MTC-0.1C) 

4  CONTINUE 
QFTUfN 
END 

RUN5E  KUTTA 

SUOSCUTIN*  KUNGE1  (T ,MECOUNT , CT , X) 

CIK^SICN  Fc(9),Pl(Q)  ,P2(9),n2<9>,X<12>,XC<12>,XT<i2), 
♦XPC12) 

T?«T 
TsTfCT 
HT*OT 
TT=T 

00  5  1*1,12 
XG (I)*X(I) 

5  XT  ( I )  *X  ( X) 

ASSIGN  6  TC 
*0  TO  2i 

6  CO  7  1*1,9 

7  XP(I>*X<I) 

«T*9.5*0T 
XTC1P)*X(IC) 

XT(11)*X(11) 

XT ( 1?) rX (12) 

ASSIGN  9  TC  K 
GO  TC  20 

9  no  1C  1*1,12 

10  XT(I)*X(I)  ‘ 

TT-T^HT 
ASSIGN  11  TC  K 

20  CALL  FOOT  <TT,PECOUNT,XT,PP) 
no  21  1*1,9 

21  X(I)*XT(I)4i.5*hT*pt(i) 

CALL  FOQT  <TT40,D»MT, ACCOUNT, X,P1) 

00  22  7=1,9 

2?  X 1 1 )  ~  y  T  (I)  C.2c7i;r7Pl*Pi(I)4.292S9??19*Pl(I)) 

CALI  FOOT  (TT*j,^*l<T  RECOUNT ,  X,  J?) 


> 


i 


"0  21  1  =  1,9 

CALL  F9CI  { TT  +  ht , T^CWNT ,  X  »  3  “ ) 
cn  3,4 

24  X(n-VTtl)  ♦HT*(c:  (I)  4»ri<*r*7S6H3S*ci(n  +  ?.i.l42125^*P2<I> 

^0  T*  «,  <6,c,ll) 

11  no  12  i*i, q 

12  X<I)=X<I)MX<I>-XP(I))/1P.1 
OFTUPN 

END 

DERIVATIVE^  FOP  Pl/JGF  KtITTA 

susrcutin*  noT  <t, -"•count,  x.xdot) 

DIMENSION  XC0T<9>  ,X<22> 

COMMON  A  A (° ,12  > 

CONSOL  L^^TICN  WCULO  "c  oPrCIFIc^  PER* 

XC1*»)  IS  ATLen',*l ,  X  (1 1 )  IS  rLEVAT0->,y<12)  IS  RUPC^ 

EXAMPLE  I  FC°  CNC  S-OCNO  DULSt  •  •  • 

IF  (T.LT.5.t,.C!?«T«Gc,6.£)  GO  10  26 
X(ll)=i.l 

no  To  ?7 

?.l  IF  (T,LT.5.c,0R.T,G“  .o.t)  no  TO  26 
X(1j)b<.,1 
GO  TO  ?7 

21  IF  (T.LT.n.O.OP.T.GF.fa.n  GO  TO  26 

V(l2)aw,  1 

GO  TO  27 

26  X(lj)aX(ll)*X(12)>I.j 

27  CONTINUE  ^ - 

no  1  J=1»9 

XOOT<J)*.).C  \K*J£2 SSSJ-5^ 

no  3  J=1 , 9 
n 0  2  K=l,12 

2  XOOT  (J)=vjCT  (J)  +A  A  (J,K)*X(K) 

3  CONTINUE 
DETU3N 
ENO 

SUBROUTINES  TO  CPTERMINF  EIGENVALUES 

SUBROUTINE  FuGVAL  <NF,  A  ,  WN,  7CTA  ,  P9,  TMALF  ,  AREAL) 
DIMENSION  A(9,6),PCR<9) ,FSI(9),7rTA  (9) ,WM(9),P0<9) 
DIMENSION  AFEAU2V 
DIMENSION  T HALF (EG 
00UpLE  PRECISION  0A{«,9> 

IF  <NC.GT.i>  GO  TO  ien 

RSKD-.J," 

GO  TO  2)C 

1<U  CO  190  T*1,N'P  W'/ 

00  190  J*1  ,NP 
19C  "A(I,J)tA(l,J) 

CALL  CHAPO  <NP,9,QAt  ..%PSR,RSI%^ 

20G  00  3:0  1*1, NP 

pomn.e  n 

WN ( I )  =  t  •  0 
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»  ^ 


TwALP(T)=j.r 
1'.*  TM!I:  ..  :  I 
:=i 

21'  IP  (I.GT.KS)  rlC  r  C 

IF  nsi(l),^,  >  GO  '*‘0  .  ,  i 

1  =  14-1 

00  TC  31? 

32f  WN (I ) =S0ST  lca-*  ( I)  **'  *-  =  51  (I)  **2) 

?ETA { T) =-RSP  f I > /WU (T)  ' 

c0lI)sf-.2i21^5,j?l^/(Wfi(I)  *S")5T  U.-7CTA 
THALF.(I)=»65V(WMI)*ZrTA(I)) 

1=1  +  2 

GO  TC  31  O' 

33:  CONTINUE 

K=1  <  ! 

00  34  0  T  =  1,N0 

TF  ANL,  °SU  (I)  •  NL  •  ‘j.  :  >  ?41t-»40 

341  A*?tAl(K)sl,/i«»S<s*3(I>)  1  • 

K=K+1  ' 

34  ■<  CONTINU* 

IF  ffcSFftLtZ)  .LT.A=>  =  U(1)  )  24  2  i  243 

342  S1=A7EAL(1)  , 

AREAL (i) =AF£AL (2)  1  1 

AorAL(2)=Sl 

343  CONTINUE  .  ’  ; 

ro  30  1=1, e 

K=I+1  1 

"0  31  J=*,c 

IF  (WN(IJ .Gt.WN(J))  PO  TO  31 
R1=WN(I) 

R2=ZETA(I) 

93=Pn { I ) 

94  =  THALcm 
WN  l I ) =WN ( J ) 

Z£TAfIJ=7ETA(J) 

PO(T)=FD(J) 

THALP(I)=THALF(J) 

HN(J)=R1 
?ET A (J) =R2 
PDtJJsR’’ 

THALF( J)=R4 
31  CONTINUE 
33  CONT+NUE 
RETURN 
FNO 

SU.OROUTINE,  CMAEC  ( N  ,  N  V  A  R  ,  A ,  C  R I T ,  c  R ,  R I ) 

THTS  SUORCUT lNr  CC‘*?UTr'S  TW£  1  EIGENVALUES  OF  A  R^Al  MATRIX, 
SYMFTRIC  C®  NONSYMMtTRJC.  Tmt  impji  MATRIX  IS  TRANSFCRMEn  oy 
FIMTIAPTTY  TRANSFORMATIONS  IN*  0  SNr  OF  ™E  FRQ^NIOUS  FORMS 
WHERE  ROW  1  CONTAIN*  ALL  OUT  THE  Lr  APING  CCEFF  OF  7HF 
CHARACTERISTIC  CCN.  T^r  LEACING  CO^FF  js  1  .  ' AC^URAC Y  TS  IN¬ 
CREASE''  0 Y  MAXIMIZING  nlVlcOR  RY  INTERCH ANdTG  "NWS  AND 
COLUMNF.TH*  ROOTS  TRC  SOLVED  USING  A  D~ AL E ^FF ^ 1  LEMMA 
WHERE  AIL  VALUES  IN  A  -"'OW  TO  THE  LEFT  CF:  TK  01  Ar-  AFC  LESS 
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THAN  INPUT  r<rjrciCot  S 'JETT  VIC*  S  P"C°  USTN'C-  °F  TO 

C'J^pfrr  n\  •>  ->?  vcrnr  L  rwr  R  POLYNOMIALS  ^M£3f  j 
A  IS  NVAk  5Y  NVAC  "ICUr3Lr  °-TC  N'  A  T  R I  X »  ,\  19  OFO'F  OF  MAI^IY, 
RR,15!  AR^  STO^AG"  AR  OA  Y  9 ,  r?jT  IS  OIVISO?  CFTT'^IA  (NC°I-'AL  G) 
CM  A  PC  USES  FOLY^F  ,LSw  =  FT,DrLYrV  AN'3  COFFER 
r'lM^9ICN  A  (‘!VAC,N»/a?)  ,  P-3  (1 )  ,  RI  (i) 

ocurle  ^cigion  nosFdi) 

DIMENSION  X  X  ( 1  ’)  ,  YYUC1 


C 

c 


3C5C 

3100 


1 

2 


,  1C 

c 

c 

c 


2C 


*0 

AT 

c 

c 

c 


A? 


DOUBLE  =::-:CIJICN  'Uw,niv,?C.<<(lLi)  ,COL  (1C) 

COL-3L c  c-rCI9ICN  A ,  X  » Yrr 

MATRIX  NORMALIZATION  FOP  A  SPECIAL  CLASS  OF  FRO?LE«<? 

IF  N  GE  20  CIV  ICE  ALL  MATRIX  ELEMENTS  py  10. C 
IF  (N.LT.21)  GO  TO  31u9 
00  3.5u  I  *  1»N 
nO  rsri  J  =  1 » N 
A  C I »  J )  =  A(I,J)/ir.C 
OONT  TNUr 
JACK  = 0 
M=N 
NR=1 
L=M 
M=L-1 
OIGsCPIT 
JJO 

FIND  LARGEST  ROW  ELEHruj  TO  LEFT  OF  DIAGONAL 
"0  lw  U=1,K 
YOD  =  A  (L  i  J ) 

AA  =  JA”S(Yr°) 

IF  (AA.Ll.SI^)  GO  to  10 

ni&  =  AA 

JJ=J 

CONTINUE 

IF  ALL  ELEMENTS  LEFT  Oc  DIAGONAL  ARE  LF  CRITERIA  GO  TO 

COMPUTE  EIGENVALUES  OF  ^EOUCEC  MATRIX 

IF  (JJ.EO.?)  GO  T0  7r 

SHIFT  FOWS  AnO  CCLS  IF  NECESSARY 

IF  (JJ.EO.K)  GO  TO  4C 

HO  2r  J=1,M 

X  =  A ( JJ, J) 

A(JJ,J)  =  A (K »  J) 

A<K,J)=X 
DO  3u  1=1, L 
X=  A ( I , J J) 

A(I,JJ)=ACIfK) 

A(I,K)=X 

CONTTNUE 

MAKE  SIMILARITY  TRANSFORMATION  CN  MATRIX 
01  V=  A  ( L  »  K) 

ROW  IN  EFFECT  ,  IS  THE  LEFT  OR  INVERSE  SIMILARITY  MATRIX 
COL  IN  EFFECT  IS  THE  RIGHT  SIMILARITY  MATRIX 

CO  42  J=1,M 
ROW  ( J)  =  A  (L  » J> 

COL (J) =-nCN C  J) /CIV 


col<m)=i. ,/rrv 
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1 


* 


0  <PCW^T)  *  A  WHFrE  KPU  IS  <TH  pnw,I  THE  IDENTITY  matrix 
DO  5C  J=i,^ 

SU^sc.CDC 
DO  45  1  =  1, F 

45  SUrt=SU.w  +  A(l,J)*FCWm 
5C  A(K,J)=«U* 

C,  <CQL*I)  *  A  WH^EC  ROW  IS  KTH  ROW, I  THE  TOENTTTY  MATRIX 

C  FIRST  K  ROWS  LESS  *TM  CCL. 

DO  6  !  1  =  1, K  . 

DO  6,  J  =  1 ,  N 
IF  (J.EO.K)  GO  TO  6? 

A(I,j)=MI,J)+MI,*£)*COL(J) 

60  CONTINUE 
C  LTH  ROW 

00  65  Jsl,” 

65  ML,J)=u.„r? 

C  KTH  COL 

A(L,<)=1. jcr 
no  63  1  =  1, K 

6*  A(I,K)=A  CI,K()*CCLCO 
L=L-1 

IF  (L.FO.l)  00  to  73 

no  to  2 

C  SET  UP  TO  COMPUTE  ROOTS  OF  RECUCEC  OR  FULL  MATRIX 
7 C  CONTINUE 

IF  (L.FO.M)  no  jn  2 J C 
COLO(l)  =  i,: 

J=1 

nO  3-  I=L,  M 
J=J*1 

COCP(J)=-A(L,I) 

3?  CONTINUE 

C  J  n-COMES  CTGPEr  OF  POLYNOMIAL 
J=J-1 

CALL  FCLVRF  (COrc,J,Xy,YY,IERR) 

TF  (IERR.NE.C)  F°INT  lOa^,  IERR 

1385  FORMAT (1HJ,1jX,13Hc0LPF  IERR  =,I3) 

C  STORE  J  PCCTS 

DO  90  1=1, J 
NR=N«H1 
RR<NR)=XX(I> 

9R  PI  (NR)  =  Y  >'{  I) 

IF  (NR.GE.K)  GO  TO  530 
M=N-NR 

IF  (M.EQ.l)  GO  to  220 
GO  to  1 

C  ONF  EIGENVALUE  IS  A  DIAGONAL  ELEMENT 

230  NP=NR*1 

PR(NR)=A(L,l)  • 
nI  (NP)  =  u.3 

21C  TF(h3«tO.N)  GO  TO  500 
TF ( L • EC . 2)  GO  TO  22^ 

KsN-NR 
GO  TO  1 
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P-wr.’  m’**F***tir#m“*  **  *** 


22* 


C 


?1 5  C 


C 

c 


5 

n 

2C 


21 


2? 

2«* 

30 


N'9=N~4j 
RR(NP)=A(t,  j) 

«?I  (NR)  *1.0 

go  to  m 

IF  N  Gc  ?.”■  MULT.  ALL  cOCTS  pY  ln.1 
IF  ( *•  .  L  T  «  2  j  )  9  "TURN 

no  jis'w  T  =  i,w 
PR(T)  =i:.,*»R(T) 

Rid)  =i:.:*rkt) 

PETUCN 

ENO 

SUBROUTINE  LEMPRT 

THIS  RCUTINf  SYSTPMATlCiLLY  FINCS  ft  RCCT  CF  ft  POLYNOMIAL 
USING  A  SIMPLE  CftGING  SCm:m£  o^Z  ON  n-ALrMnERTS  LEMMA 
COMMON  /COFFES/F^dl) ,  M,  X  ,Y  ,  ftp ,  ®p,  T I  ,IE9  ?R 
01  PENSION  NPLAG(C)  ,0  (5)  ,  V  (5  )  ,p(5> 

DOUBLE  PPEC  T5I CN  P* 

EQUIVALENCE  (F  , F 1 ) ,  ( F  ( 2 ) , F2) ,  (F  (?) , F ?> , (P (4 ) , P4) , 

*  (P  ( 5) , P5) 

L  =  1 
9R  =  :.** 

pi  * 

SIGN  a  1.0 
IFLAO  *  0 
JFLAG  =  0 
KFLAG  s  Cl 
PEL  *  C.? 

POLL  «  «.0 
CP  51=  1,5 
MFLAG(I) s  0 

IF  (IcL4G.LT,5)  GO  TP  25 
IF  (JFLAG.GT.S)  GO  TO  ?5 
IFLftG  s  « 

IF  (YFLftG.LT. 3)  GO  TO  22 
RR  =  9P  +  STGN/1® .  C 
91  a  RI  +  «IGN/1?.0 
*IGN  =  -2.‘J*9IGN 
CONTINUE 

IF  (flPS(PI)  .LE  ,1.0  GO  to  22 
SIGN  =  SIGN/97.0 
R9  s  STGN/3.0 
RI  =  -SIGN 
CO  TP  21 

KFLftG  =  KFLAG  4  1 
DEL  *  nOEL*CEL 
OOFL  =  OOEL  4  1.3 
►.'FLAG  (L )  =  C 
GO  TP  ?; 

IFLAG  =  IFLAG  4  1 
CONTINUE 
CO  41  I  *  1,5 

IF  (NFLAG(I).NE.C)  CO  TO  38 
X  =  RR 
Y  =91 
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!F  C I  •  L  ^  *  I )  go  n.  *»«? 

IF  X  r  x  f  [ti 


IF  C.tO,?)  X  r  V  .  ^cL 

IF  (I.^.L)  Y  *  V  ♦  CTL 

IF  (I.LO.F)  Y  =  Y  -  "EL 

•*5  U(I)  =  X 
V  <  I  >  =  Y 
CALL  FCLYEV 
cm  r  as 
NFL  A" (I)  s  r 
4o  CONTINUE 

IF  (JFLAG.GT.27)  GC  TO  6* 

°0  «'  I  ■  ^ 

IF  (P(I)  ,GT «  1.;*-j7>  GO  TO 
4  5  CONTINUE 
go  y c  f ; 

4e  0IF1  =  AMAXl(Fl,c2»P’i,fp4»P5>  ^ 

PIF2  =  AMIM(Pl,C2,P?,p4,P5) 

DIP  v  Clpl  -  0IF2 

IF  ((  CIF.Gf.i.C).4NP.(Pl.LT.l*'}>)  GO  to  *«? 
IF  (=i.E0.u.->  GC  TO  60 
DIF  s  PTF/ci 

IF  OIF. LT. 1.001)60  2T 

55  CONTINUE 
6r  rONTTNUc 

CO  7:  J  S  1,5 
I  s  J 

IF  (P(J).fc.-.w)  GO  TO  1.0 
70  CONTINUE 

"IF2  ■  AYIM(P?,P?,P4,P5) 

IF  (F1.GT.CIF2)  GO  TC  8C 
IF  (DEL.LT.n,  j£-33)  RETURN 
PEL  =  C.5*CEL 
XX  s  pk  ♦  CEL 
YY  s  RI  ♦  CEL 

IF  (  (XX.FQ.p*)  .AND.  (YYtEF**^!) )  RETURN 
IF  ((XX. EO.P3) .AND. (PI. LO. ..:))  RETURN 
IF  ((Ro.EO.C.d.AMC.tRI.FQ.YY))  RETURN 
IF  (JFLAG.GT.l'w)  GO  TO  220 
JFLAT,  =  JFLAG  ♦  1 
NFL  AG  (1)  =  1 
GO  TC  2j 
3C  AMINY  *  02 
N  =  2 


DO  35  1  =  7,5 

IF  <FU)  .GT. AMINY)  GO  TO  85 
N  *  I 

AMINY  =  P(I> 

*5  CONTINUE 
L  =  ? 

IF  (*.50.3)  L  -  2 

TF  (M.10.4)  L  a  5 

IF  L  *■  4 

NFL  AG (J)  =  1 
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NFUG(L)  =  i 
UCL)  =  U(l) 

U  ( i )  =  M(M) 
VCL)  =  V ( i  ) 
P(L)  =  PI 
VC  1  >  .=  V  (N) 
PI  =  P(M) 
RP  =  UC1) 
rI  *  V  ( 1 ) 

r’0  to  13 

1-C  RR  a  UCI) 

PI  =  VCI) 

return 

22 ii  IERRR*  2 


RETURN 

eNO 

SUBROUTINE  FOLYEV 


COVMCN  /CC  =  rER/P='(ll)tM,x,Y,A.F,Ri?,RI 
OOUTLE  °Rp,CI5TON  PR,U,V»US 
U  «  fiR  ( 1 ) 

v  *  j. "on 

DO  2.  I  *2,K 
US  =  U 

U  =  X*U  -Y*V  ♦  P  R  ( I ) ": 

2C  V  =  X*V  ♦  Y+U3 

AP  a  OARS(U)  ♦  C5PSCV) 

RETURN 

EM? 

SUP  ROUTINE  FOLYRF  (°» M, X , Y , IEPR) 

CC^'FCN  /COcF~R/CR  til)  ,M,A,E>ftF,RP,RI,IFRS{F 
00U9LE  PRECISION  PCI)  ’ 


RIMlNSICN  X  (1) ,YC1) 

D0U9U  PRECISION  PR,P,X?,XY 

IF  t (N.LT.l) .OR, (N.GT.RO) )  GO  TO  2J« 

IERR  a  0 


IERRR*  I 


J  e  1 


M  s  KJ  +  1 
00  lg  !  8  1 , M 
PR (I)  *  P(T> 

1C  CONTINUE 

IF  CN.FQ.l)  GO  tq  l'„r, 

IF  CONTINUE 
CALL  LEMPRT 

IF  (RI.EO.C.j)  GO  TO  40 
IF  (’?,E0,c.3)  GO  TO  16 
TEST  a  RI/RP 

IK  (APSCTEST)  .LT.?. 3*00*1)  GO  TO  43 
16  CONTINUE 

XCJ)  *  RR 
YCJ)  *  RI 

IF  (IfRRRDT.Q)  *<0  TO  220 
IF(J.EO.N)  RETURN 
J  =  J  ■»  1 
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23  0 
1200 


XU)  *  -3 
Y(J)  =  -31 
TF(J.?r.V)  trx y*K 
J  =  J  ♦  1 
M  S  V  m  2 
X2  =  2 .  ^  »S=» 

YY  S  -(»R*P-->  + 

<1(1  21  T  s  f,y 

FR(I)  e  Cj(j)  +  x?»P^fT-l) 

cR<I+i)  s  pMI  +  l)  ♦  YY*PR(T-1) 

20  CONTINUE 

IF  (M.LQ.?>  GO  TO  10 i 
GO  Tf  1« 

A r*  CONTINUE 
X(J)  =  RS» 

V(J)  r  3.3 

IF  (IFK^P.Nf,^)  no  TO  2?C 
IF  (J.LT.N)  RFTU3N 
J  8  J  ♦  1 
M  s  M  -  1 
CO  *0  I  s  ?,M 

5t  PR ( I )  s  Fk ( I )  +  c F+PF ( I-l ) 

IF  <M.C0.2)  no  TO  ICC 
G0  Tr  15  ‘ 

5t  0  *  CR(l> 

x(j>  =  -pp(?)/n 

V(J)  a  J.Q 
RETURN 

50  PRINT  12U,  N 

jC  FORMaTaHo,l"Xf  3HN  s,  14,  RiHOUTSIRE  LIMITS  POLYRF) 


IFRRsl 

RETURN 

220  TERR  s  IERRR 
RETURN 
FND 

DATA  CARO  PO?  CONSTANT  TERMS 
32.174  1745.5  132,6  13.7 

END 


13.71  3.1416  1  <1 2  C  •  0  10.  P9 
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a  degree  of  Bachelor  of  Science  and  a  commission  in  the 
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